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We demonstrate how most of the light J p = + spectrum below 2.0 GeV and their decays can be 
consistently described by the unitarized quark model incorporating the chiral constraints of Adler 
zeros and taking SU(3) breaking effects into account. These resonances appear as poles in the 
complex s plane in a unified picture as qq states strongly dressed by hadron loops. Through the large 
N c analysis, these resonances are found to naturally separate into two kinds: o, k, /o (980), ao (980) 
are dynamically generated and run away from the real axis as N c increases, while the others move 
towards the qq seeds. In this picture, the line shape of ao(980) is produced by a broad pole below 
the KK threshold, and exhibits characteristics similar to the a and k. 
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I. INTRODUCTION 

The enigmatic spectrum of light J p = + scalar res- 
onances are of great interest for its importance in un- 
derstanding chiral symmetry breaking and confinement 
in QCD. Despite many theoretical efforts, the current 
understanding of the microscopic structures of these res- 
onances is in a well-known unclear situation as summa- 
rized in Particle Data Group (PDG) [lj: qq models 0, 
the unitarized meson model Q , a tetraquark model with 
and without qq mixing [J, |5| , the Julich meson exchange 
model @ , the unitarized a model [7| , ghieball (8j or using 
the inverse amplitude method (IAM) ||, NJL model pjf 
and lattice simulations pT|. and so on. Most of these 
studies focus on the lowest putative nonet or explain the 
lighter and heavier resonances in different ways. In the 
present paper, we show that, all the light scalar spec- 
trum below 2.0 GeV except a glueball candidate can be 
described (or even predicted) using just seven parame- 
ters in a unified and consistent picture, that is, qq seeds 
strongly dressed by hadron loops. The picture brings 
more insights on the origin of the resonances, which are 
generated as the poles of the S matrix and have no one- 
to-one correspondence with the nonet in the Lagrangian. 
At the weak coupling limit as N c increases, <r, k, ao(980), 
and /o(980) move away from the real axis on the com- 
plex energy plane, whereas all the other heavier J p = + 
resonances move to the bare seeds. This reveals the dif- 
ferences between the lighter mesons and heavier states. 

We use the unitarized quark model (UQM) [2| pro- 
posed by Tornqvist, which played a pioneering role in 
the resurrection of the a meson. The merit of this model 
is that it naturally respects the unitarity of the S-matrix 



but also incorporates some dynamics at the same time. 
Besides, the Adler zeros [12J, as the constraints from chi- 
ral symmetry, can also be easily implemented. Neverthe- 
less, the k resonance was not found in his explicit analysis 
of experimental data, and those resonances with higher 
masses than 1.5 GeV are also not covered. In the present 
paper, however, by incorporating the SU(3) breaking ef- 
fects in the coupling constants and analyzing the poles on 
the complex plane, we show that the k resonance can re- 
ally be found in this picture. Moreover, most resonances 
in I = 0, / = 1/2, and 1 = 1 channels can find their 
corresponding poles on the complex plane. 

The paper is organized as follows: In Sec. II, we briefly 
introduce the basic scheme of UQM and the three non- 
trivial improvements we make to this model. Our numer- 
ical results are elaborately discussed in Sec. III. Section 
IV is devoted to a further study on the characteristics of 
these resonances based on the large N c technique. Sec- 
tion V summarizes our main results. 



II. THE THEORETICAL SCHEME 

A. Unitarized quark model 

The unitarized quark model begins by assuming that 
there are qq bare bound states generated in QCD and 
they are coupled with the pseudoscalar mesons. The 
main idea is to take into account the hadron loop dress- 
ing effect in the propagators of the bare qq states @, LL3| ■ 
The bare propagator of a qq bound state is 



P = 



1 



(1) 
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where too is the bare mass. For example, too is 2to for 
uu or dd, rh + Am for us, and 2to + 2Am for ss state, 
respectively. The vacuum polarization function, II(s), 



which represents all the possible two pseudoscalar meson 
loops, will contribute to the full propagator as 



P = 



1 



m.Q — s + n(s) ' 



(2) 



As an analytic function with a right-hand cut, its real 
and imaginary parts are related by a dispersive integral 



1 r°° 

ReII(s) = -V / dzImn(z)/0 - s), 

T Js th 



(3) 



where 



Imn( S ) = -]TG 4 ( S ) 2 

i 

= -Y j g} k ^-F l {sf6{ S - Sth ^ (4) 

i 

where the general coupling function Gi(s) includes the 
coupling constants gj's, the phase space factor ki{s)/y/s, 
and a Gaussian form factor Fi(s) = eocp[— kf(s)/2ko\. 
fej(s) is the ith channel cm. momentum with fcj(s) = 



A(s,m^.,m^.)/4s and the 9(s — Sth,i) is a unit step 
function. 

If there exists more than one bare state in the i — > j 
channel, the partial-wave amplitude can be represented 
in a more general matrix form: 

{P~ 1 } a /3(s) = (m^ a - s)S a p + n a/ s(s), 
Imn^( S ) = -VG M ( S )G^( S ) 



V^ " ■"» i ' s ^ t? 2 



g a i9pi jL i=-F i (s)6(s - s t h,i), (5) 



where ReII a/ 3 is determined by a similar dispersion inte- 
gral of Imilo,^ as Eq. © . The off-diagonal terms of H a p 
produce the mixing between different bare states coupled 
with the same intermediate states. 

The Adler zeros are incorporated into the UQM model 
in a direct and easily operated phenomenological way [2j : 

k(s) 
G a i(s)Gf3i(s) -> 7 Q i7/3i(s - z A}i )F 2 (s)- ! -^9(s - s thA ), 

(6) 

where the z A /s denote the Adler zeros, and 7^ are di- 
mensionless coupling constants . 



B. The Adler zeros in the Chiral Perturbation 
Theory 

Normally, the T matrix also contains left-hand cuts. 
Because the Adler zero is usually located nearer to the 
physical threshold than the left-hand cut, it is natural 



to expect that it plays a more important role than the 
left-hand cut in determining the scattering amplitudes 
along the right hand cut. Since such zeros reflect the 
constraints of the chiral symmetry, in the / = 0, 1/2 
channels, they are fixed in our study at the values ob- 
tained from the chiral perturbation theory (ChPT) [14( 
rather than left as free parameters. In the 1=1 chan- 
nel, we also fix the Adler zero according to ChPT but 
there are some subtleties which will be addressed later. 
Being within the convergence radius of the chiral expan- 
sion, these Adler zeros should be reliably determined by 
ChPT and, hence, is a reasonable starting point for a 
phenomenological study. 

After partial-wave expansion, one obtains the Adler 
zero of I = tttt S wave at about rn^/2 at tree level, and 
the position is slightly shifted to s ~ 0.38to 2 by including 
the contribution up to two-loop SU(2) ChPT15]. Sim- 
ilarly, in the Kn scattering, the SU(3) ChPT to 0{p 2 ) 
gives S-wave amplitude of I = 1/2 as 



rjz? /a (*) = 



-3(m K 2 



2 {rriK 2 + Tow 2 ) S + 5s 2 



128nf K ^ S 



(J) 



which has two zeros located at ^(m,K 2 + rn^ 2 ± 
2\fArriK i — Ittik 2 fn-K 2 + 4rn w 4 ) in the unphysical region. 
One is on the negative real axis inside the circular cut, 
and the other is on the positive axis between the circular 
cut and the Kn threshold. With the 0{p A ) contribution, 
the left zero on the negative axis moves to the complex s 
plane at about s = 0.003 ±0.005zGeV 2 and the right one 
on the positive axis will also be slightly shifted (at about 
s = 0.233GeV 2 ). As for the ixr\ scattering, at leadin g or - 
der 0(p 2 ), ChPT recovers the current-algebra result [la]: 



T$(s,t,u) 



3/ 2: 



(8) 



which contains no zero on the real s axis after the partial- 
wave projection. Including the higher order contribution 
is not helpful to obtain a real- valued Adler zero. A pair of 
zeros, at about 0.078 ± i0.178 GeV 2 , can be found when 
the 0(p ) terms are taken into account using the low- 
energy constants from Rcf.[17]. Since there is no accurate 
data of the irr] scattering, we could choose a zero point 
on the real axis to simulate them as mentioned later. 



C. The scalar- pseudoscalar-pseudoscalar coupling 

The coupling of pseudoscalar-pseudoscalar to the + 
states could be described by effective interaction terms in 
the Lagrangian: L S pp = aTr[SPP] + /3Tr[S]Tr[PP] + 
^Tr[S]Tr[P)Tr[P}. The first term has an SU(3)- 
symmetric quarkonium coupling as used in many phe- 
nomenological models. We also include the Okubo- 
Zweig-Iizuka violation terms in the last two terms in a 
general way. Moreover, the decay of the quarkonium into 



TABLE I. The effective scalar quarkonium coupling to pseu- 
doscalar mesons up to a global constant. 
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FIG. 1. The right-hand cuts and their signatures in a three- 
channel case. 



D. The analyticity of S matrix 

By definition, a resonance is specified as a pole of the 
S matrix analytically continued to the complex s plane. 
Extracting the poles of the partial-wave amplitude of 
i — >• j process described in the UQM is actually to find 
the zeros of the determinant of the inverse propagator. 
Sometimes, it can be obtained in some other equivalent 
way. For example, in a two-channel occasion the analyt- 
ically continued S matrices on different Riemann sheets 
could be written down using those on the first sheet [19j: 



a pair of mesons QQ — > M(Qq~i)M(qiQ) involves the cre- 
ation of a qiq~i pair from the vacuum. The ratio of the 
creation rates of ss and uu or dd from the vacuum is usu- 
ally defined as p = (0|V| ss) / (Q]V\uu) , representing the 
breaking of SU(3) symmetry [IS]. SU(3) breaking ef- 
fects have proved to be important, so we allow for these 
in our version of the UQM. To be explicit about our de- 
scription of the coupling of quarkonium to mesons, we 
express the scalar and pseudoscalar 3x3 flavor matrices 
in qq configurations as 



^a° + -fef n a+ n+ 



K 



V2 J V2 J 



fs 



/ V3ir +y 8 



P 



"K 

K 



7T+ K+ 

-V?,-K a +r l& ts-0 

V6 



K u 



hs J 



(9) 



^,(10) 



where /„ = nn = (uu + dd)/\/2 and f s = ss. The 
physical states, r\ and 77', are conventionally defined as 
X] = coscf)\nn) — sin0|ss), 77' = smc/)\nn) + cos</>|ss), with 
cf) = tan _1 \/2 + Op, where dp = —11.5° being the pseu- 
doscalar octet-singlet mixing angle [l|. Thus, by stan- 
dard derivation, a general form of effective coupling con- 
stants between the scalar quarkonium and pseudoscalar 
pair is obtained, as shown in Table HI 



S n = 



?iv 
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iS 12 \ 


Su 


Su 


iS 12 


detS 


Su 


Su / 


detS 


iSi 


S22 


S 2 


1S12 
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S22 


S22 



,s 



III 



S22 

ti< IS 


—S12 
detS 


— S12 
detS 


Su 

del S 



(11) 



which implies that a pole on the second Riemann sheet 
is just located at a zero point of Su on the first sheet, 
a third-sheet pole at a zero point of detS, and a fourth- 
sheet pole at a zero of S22, respectively. 

In the literature it is common to define a Brcit-Wigner 
mass of a resonance as the solution of rn BW = m B + 
ReH(m BW ). This is a good approximation for narrow 
resonances and is commonly used in experimental anal- 
ysis. However, if the propagator is strongly dressed by 
hadron loops where 7mII(s) is large, the mass and width 
are no longer suitably determined by the Breit-Wigner 
form, but should only be defined by the pole position of 

Sp + n(s p ) = 



'B 



the S matrix, i.e., the solution of ■ 

with s p = (M — iT/2) 2 . Because of the analyticity of the 
S matrix, the determinant of inverse propagator vanishes 
only on unphysical Riemann sheets. 

The general character of the poles on different Rie- 
mann sheets has been discussed widely in the literature, 
(see, for example, [2fJ,[21(). Every physical cut will dou- 
ble the Riemann sheets in the analytical continuation, so 
there are 2™ Riemann sheets in a process with n coupled 
channels, as shown in FigfT] The physical sheet is defined 
as the sheet where all the cm. momenta are positive on 
the physical cuts, denoted as (+ + + •••+) signature. In 
the same fashion, the (n + l) th sheet (n < N), attached 
to the physical sheet between $th,n and Sth,n+i along the 
physical cut, is denoted by (— • • ■ — !- + •••) with the n 



TABLE II. The results of the fit parameters with a Gaussian 
form factor, mo, mo + "ij, and mo + 2m s are the bare masses 
of nn, ns, and ss respectively. 



a : : 7 


1.493 ±0 .05i : (-0.149±o.02 5 ) : 0.319 ±o .02i 


P 


0.704 ±0 .054 


fco(GeV) 


0.505±o.oo9 


mo(GeV) 


1.443±o.o2o 


m s (GeV) 


0.046±o.oi2 



consecutive "— " signs before the other "+" signs. A res- 
onance is represented by a pair of conjugate poles on the 
Riemann sheet, as required by the real analyticity. The 
micro-causality tells us the first Riemann sheet is free of 
complex- valued poles, and the resonances are represented 
by those poles on unphysical sheets. The resonance be- 
havior is only significantly influenced by those nearby 
poles, and that is why only those closest poles to the ex- 
periment region could be extracted from the experiment 
data in a phenomenological study. Those poles on the 
other n > N sheets, which are reached indirectly, make 
less contribution and are thus harder to determine. 
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FIG. 2. The left three figures show the fit quality. The right 
three ones are predictions. Compared data are from 25, 26]. 



III. NUMERICAL ANALYSIS 

Now, we apply the partial-wave formulation 2] with 
our new ingredients to study / = 1/2 Kir, I — 7T7t, 
and I = 1 irrj S-wave scattering. The main purpose of 
this paper is not to make an exhaustive fit, so the only 
data used in the combined fit are: (1) the I = 1/2 S- 
wave Ktt scattering amplitude I22ll23ll . (2) the phase shift 
of J = S-wave 7T7t scattering [241 ] . and (3) the phase 
shift of tttt -► KK [H below 1.5 GeV. The fit quality 
is good and x 2 /d.o.f. ~ 0.8. The central values and 
the statistical errors of the seven parameters are listed in 
Table |TTJ The good agreement between our theoretical 
results and experimental data can be seen in Fig. [5J even 
though some of these data have not been used in the 
fit. The parameter values in Table [IT] are all in realistic 
ranges. The bare masses of qq states are slightly larger, 
but not in conflict with NJL modelings [13 ■ The SU(3) 
breaking effect parameter is also consistent with the value 
in the literature [18.]. A general comparison of the masses 
and widths of the resonances from our results and the 
values from the PDG table is presented in Fig. [3] which 
will be discussed in detail. 

The 7 = 1/2 Ktt S-wave scattering provides an ideal 
illustration and the best testing ground because of its 
large threshold spans and clean experimental informa- 
tion. With the parameters of Table [TT1 there is only one 
solution to s = (too + m s ) 2 + ReII(s) for real s at about 
s 1 / 2 = 1.33 GeV, which is the Breit-Wigner-like mass 
mentioned previously. However, there are three poles of 
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FIG. 3. The filled symbols represent the resonances' masses 
and widths from the PDG table, and the empty ones represent 
the pole masses and widths we obtained. The reason of the 
discrepancy between the value sets has been discussed in the 
text. 



the S matrix found near the physical region at 

V7 1 = 0.767 ±0 .oo9 - ^0.308 ±0 .o35, 

V7 77 = 1.456±o.ois - z0.164 ±0 .o26, 

V?v = 1.890 ±0 .o29 - z0.296 ±0 .oi4, (12) 

where the superscript denotes the number of the sheet 
and the units are in GeV. Simply comparing these poles 
with the tables of Particle Data Group [l| , a good agree- 
ment in quality is instantly found (see Fig. \3§ . The lowest 



second-sheet pole is just the n resonance and consistent 
with the values determined by those mo del- independent 
methods [23]. The third-sheet pole corresponds to the 
^(1430) and a second-sheet "shadow" pole (due to the 
weak coupling constant of the Krj channel in our result 
and also found in [21]) is also found at almost the same 
location. Although the fit is only carried out with the 
data below 1.5 GeV, a fourth-sheet pole is predicted by 
the SU(3) couplings and unitarity constraints, which cor- 
responds to the higher Kq(1950) resonance. The width 
of iCg(1950) is larger than its PDG value, but is quali- 
tatively acceptable compared with the average value cal- 
culated from both solutions A and B of the original data 
analysis pl[28J|. 

The poles in / = 7777 S wave are inevitably more 
complicated than those in the Kn S wave because of the 
mixing of nn and ss states. We find 

V7 1 = 0.430 ±0 .o4o - i0.249 ±0 .o75, 
~U _ 

a™ = 1.467 ±0 .035 
/W = 1.577 ±0 .o4o- 

VI 



TABLE III. The results of the fit parameters with another 
form factor. 



's" J = 1.935 ±0 .028 
W 1 = 2.444 ±0 . 03 2 



i0.228±o.064, 
i0.306 ±0 .o23, 
i0.289 ±0 .oi3, 
i0.242 ±0 . 013 . 



(13) 



All these poles could be assigned to those light resonances 
of I G (J PC ) = 0+(0++) listed in the PDG table, except 
for the /o(1710). The position of the a pole is in agree- 
ment with the results of model-independent analysis [29| , 
while the /o(980) is a narrow pole below the KK thresh- 
old. The /o(1370) is a fourth-sheet pole and its position is 
within the uncertainty of the PDG value. The pole mass 
is consistent with that preferred by the Belle Collabora- 
tion from 77 — > 7r°7r°, at 1.47GeV [3fJ. Here, the res- 
onance shape of /o(1500) is generated by the fifth-sheet 
pole and 7777' threshold together, as found in Ref. [3lj | . 
This may be the reason why the width of the pole is 
much wider than the PDG value. The Air and pp thresh- 
olds turn out to be increasingly important beyond about 
1.2 GeV. However, the inclusion of such thresholds re- 
quires the SVV and VPP interactions to be taken into 
account. This would introduce many new parameters and 
this case is beyond the scope of this paper. Not incor- 
porating the SVV interaction might explain why there 
is no Jo (1710) pole in this picture. The other possibil- 
ity is that the main ingredient of /o(1710) could be the 
lowest scalar glueball, as preferred by recent quenched 
Lattice calculations [32j. This would add a narrow res- 
onance structure with its own hadron cloud [33]. The 
two sixth-sheet poles are assigned to the /o(2020) and 
the /o(2330), respectively. 

As for the I = 1 ~ki) scattering, owing to the poorness 
of data, this channel is not included in the fit and so 
the plot in the lower right corner of Fig. 1 is wholly a 
prediction. As we have mentioned previously, no real- 
value Adler zero is found up to 0(p A ) ChPT amplitude 
in this partial wave. We simply set za = 0.078 GeV 2 in 



Q : /3 : 7 


1.889± .o24 : (-0.280 ± ,,.oi5) : 0.269 ±0 .on 


P 


0.821±o.O27 


mo (GeV) 


1.763±o.o38 


m s (GeV) 


0.215 ±0 .037 



the calculation, which is close to the complex zeros. It is 
not difficult to exhibit the ao(980) line shape below the 
KK threshold, as shown in Fig. [51 The pole positions are 
not sensitive to small deviations from the value of Adler 
zero we choose: 



Vs 11 = 0.792 ±0 .oi5 - z0.292 ±0 .o6o, 

Vs 171 = 1.491 ±0 .034 - l0.133 ± o.O38, 



1.831± .o27-*0.265± 



0.014- 



(14) 



The second-sheet pole below the KK threshold is broad, 
but still produces the ao(980) line shape combined with 
the threshold effect, as proposed by Flatte [34[. This 
effect was also found by studying the 7777 amplitude in 
the Julich model [351 ] and implied in the unitarized a 
model [7J. The third-sheet pole could represent the 
ao(1450), although it plays a minor role in our picture. 
There is also a higher ao(1830) predicted, which has not 
been widely observed in experiments, but might be re- 
lated to the ao(2020) seen by the Crystal Barrel Collab- 
oration [361 ]. 

The use of a Gaussian form factor is the most dras- 
tic assumption we have made. This has been widely 
used in experimental analyses and in many other mod- 
els. However, an exact representation of the form factor, 
satisfying unitarity and analyticity and easily applied in 
phenomenological studies, has not been found. To test 
the stability of our results, we also use a different form 
factor, ](M 2 + s th ,i)l{M 2 + s)] 2 (M is the mass of the 
resonance) proposed in [33]. The latter form factor is 
not better than the Gaussian form factor since it suffers 
from a spacelike pole, but it might provide a reasonable 
qualitative cross-check, especially for the N c analysis we 
addressed later. Incorporating the latter form factor re- 
moves one parameter, ko, so the fit quality is worse and 
X 2 /d.o.f. ~ 1.21. The central values of the other six 
parameters listed in Table IIIII are different from those 
in Table llll Nevertheless, the poles below 1.5 GeV, as 
shown in the following, coincide in position with those 
found using the Gaussian form factor: 



0.745± .oo7-«0.301 ±i 



0.023, 

V7 77 = 1.547 ±0 .o2i - 70.148 ±0 .o26, 
for the / = 1/2 channel, 

VI 77 = 0.396 ±0 .oi9 - i0.244 ±0 .o46, 
V7 1 = 0.984 ±0 .o23 - i0.028 ±0 .o26, 



(15) 



for the 1 = channel and 



0.700 



±0.060 



i0.265 



±0.031, 



(16) 



(17) 



for the 1=1 channel. The corresponding poles beyond 
1.5 GeV still exist, but since we only fit the data below 
1.5 GeV and have not included the 1=1 data, they move 
farther away from the previous values. A posteriori, this 
means the second kind of form factor may not be a good 
choice. 
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FIG. 4. Left: 1=1/2 poles' trajectories; Right: the pole tra- 
jectories of a and /o(980). 



IV. LARGE N c ANALYSIS OF THE POLE 
TRAJECTORIES 

The success of describing such a broad range of spec- 
trum and their decays in a unified and consistent way 
suggests that it is a reasonable model to study these res- 
onances and could be used to gain further insights into 
their nature. The large N c behavior of the pole trajecto- 
ries serves to shed light on the origin of these resonances. 
The lowest order of a is l/\/W~ c , (3, and 7 by a factor of 
1/N C and 1/N^. The bare mass, the location of Adler 
zero, and the form factor are of order 1, while the mixing 
angle <fi is of order 1/N C [3Ej. Whichever of our two form 
factors we use, the poles exhibit similar trajectories as 
iV c increases. Those for the I = 1/2, resonances are 
shown in Fig. @] as examples. The cr, k, and ao(980) 
poles move farther away from the real axis. In contrast, 
the Kq(1430) and A~q(1950) become narrower and move 
towards the ns bound state. Analogously, the ao(1450) 
and ao(1830) move to the (uu — dd)/V2 bare seeds. The 
/o poles other than a and /o(980) move towards either 
the nn or ss bare seeds. At N c = 3, if the coupling 
to the 7T7T channel is switched off, the /o(980) will move 
down below the KK threshold and form a bound state. 
This behavior implies that the /o(980) is more like a KK 
molecule state. However, when N c increases, it exhibits 
a peculiar trajectory: it moves to the real axis rapidly 

and then crosses the cut onto the (H h + +) sheet, and 

then moves away from the real axis as the a pole does, 
as seen in Figf4j If the coupling to the lowest thresholds 
are increased, respectively, by hand when N c = 3, the a, 
k, and ao (980) will move to the real axis and become vir- 
tual bound states different from the seeds either. While 
the coupling becomes strong enough, the virtual bound 
states will move onto the first sheet and become bound 
states. It is worth mentioning that, in using the IAM 
to unitarize ChPT [39] , Pelaez has observed similar pole 
behaviors of a and n in some parameter region. Using 
the Pade technique to unitarize ChPT amplitudes, the 
similar a pole trajectory is also found 40] . The pole be- 
haviors of Jo (980) and ao(980), as we pointed out here 
for the first time, may explain the strange behavior of 
the line shape in large N c shown in [39j. 

So, the general N c behavior separates the poles into 
two types: a, n, ao(980), and Jo (980) are the first type (or 



the unconventional type) of resonances, like bound states 
of mesons, which are dynamically generated by the pseu- 
doscalar interactions. This may be the reason why they 
could be described by the tetraquark model. All the other 
resonances except the glueball candidate, as the second 
type (or the conventional type) of resonances, are directly 
generated from qq seeds by renormalization effect, which 
indicates that they all belong to the same bare qq nonet. 
As the interactions are turned on and different channels 
are open, the bare seeds are copied to different Riemann 
sheets and get renormalized by the hadron clouds in var- 
ious ways [2l| . Some of them run too far away from the 
physical region to be detectable. In this picture, there is 
no need to distinguish parts of them to be a nonet. 



V. SUMMARY 

In conclusion, this paper demonstrates that the whole 
low-energy scalar spectrum below 2.0 GeV, except for 
a possible glueball /o(1710), could be described in one 
consistent picture, with the bare u qq seeds" dressed by 
the hadron loops. All the resonances are dynamically 
generated by the same mechanism, and there is no di- 
rect correspondence between the poles and the original 
nonet in the Lagrangian. In a large N c analysis of this 
picture, the pole trajectories exhibit a general behavior 
which agrees with other models. In particular, the c, k, 
/o(980), ao(980) resonances, though running away from 
the real axis when N c is larger, are also generated in this 
model, which means this large N c behavior does not con- 
flict with the qq dressed by the hadron loop picture. They 
are produced by large hadron loop effects and this may 
also imply their large four-quark components. Thus, in 
this paper, we present that the usual speculation in par- 
ticle physicist community, that the lighter scalars behave 
like the tetraquark states and the heavier scalars do as 
the qq states, could be actually realized in such a coherent 
picture of improved UQM model. 

We also show how the line shape of ao(980) is possibly 
generated by a deep pole, like the a or k, encountering the 
KK threshold. This whole treatment could be extended 
to other spectra, e.g., the charmonium states [41], and 
provide theoretical suggestions for further experimental 
investigation. 
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